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Abstract 

The shear viscosity of a two-dimensional liquid-state dusty plasma was measured experimentally. 
A monolayer of highly charged polymer microspheres, with a Yukawa interaction, was suspended 
in a plasma sheath. Two counter-propagating Ar + laser beams pushed the particles, causing shear- 
induced melting of the monolayer and a shear flow in a planar Couette configuration. By fitting 
the particle velocity profiles in the shear flow to a Navier-Stokes model, the kinematic viscosity was 
calculated; it was of order 1 mm 2 s , depending on the monolayer's parameters and shear stress 
applied. 
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A broad range of charged particle systems can be modeled as Yukawa systems that 
can have a liquid state; these include colloids, certain dense astrophysical plasmas, and 
strongly coupled dusty plasmas. A dusty plasma is a suspension of highly charged micron- 
size particles in a plasma. When these particles are confined, their mutual repulsion causes 
them to self-organize in a structure called a plasma crystal, which can be in a crystalline or 
liquid state. 

Shear viscosity is a dynamic property of fluids required to describe shear flows and damp- 
ing of waves. Recently, molecular dynamics (MD) simulations and theory [4] have 
been used to predict the viscosity of a liquid Yukawa system. While viscosity is commonly 
measured in colloids, there are no reports of comparable measurements in dusty plasmas. A 
dusty plasma is much softer than a colloid, and between its particles it has only a rarefied 
medium consisting of free electrons and ions, as well as neutral gas, which applies a frictional 
drag to moving particles. 

A plasma crystal is an analogue for molecular matter; the particles represent molecules 
0, but they have the advantage of allowing direct imaging and thns a n_nent of 
their individual positions and velocities pj. Confinement of particles is provided by natural 
electric fields in the plasma. The particle suspension can be manipulated using various forces. 
In Ref. ion drag forces were used to apply a shear stress that melted a two-dimensional 
(2D) crystallized dusty plasma. In Ref. p], a 3D dusty plasma in a liquid state flowed 
past an obstacle where there was a shear. In those two experiments, experimenters had no 
direct control over the forces and the resulting shear. Here, we report an experiment using 
laser radiation to apply a controlled and localized shear stress to a monolayer lattice. The 
lattice melted, and we modeled the resulting velocity profile to measure the nonequilibrium 
viscosity. The geometry resembled a laminar planar Couette flow. 

We used the apparatus of Ref. 2], with essentially the same parameters, including an 
Ar pressure of 5 mTorr. A monolayer of microspheres was suspended in the plasma. The 
particles had a diameter of 8.09 ± 0.18 /iin jj| and a mass m = 4.2 x 10~ 13 kg. The 
particle suspension's diameter was 50 — 60 mm. The interparticle potential for particles 
arranged in a single plane, like ours, was experimentally shown to be nearly Yukawa 
U(r) = Q(47re r')~ 1 exp(— r/\ D ), where Q is the particle charge and is the screening 
length. (This should not be confused with quasi-2D arrangements of particle chains, which 
might have a different potential.) The particle suspension is characterized by k = a/Xa and 
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r = Q 2 /4ireoakT, where T is the particle kinetic temperature. For liquids, the characteristic 
length a is the 2D Wigner-Seitz radius [11]; it is related to the lattice constant b for a perfect 
triangular lattice by a = (y/3/2n) 1 ^ 2 b. To vary a, k, and Q, while keeping the plasma 
parameters constant, we used a different number of particles in each of three experiments, 
Table 1. We used the pulse technique of Ref. jl2| to measure k and Q. 

Initially, our suspension was an undisturbed triangular lattice, Fig. [TJ^a). As is typical for 
a 2D monolayer with our experimental conditions, it had a highly ordered state. The static 
structure factor S(k) has the distinctive peaks of a hexagonal crystal, Fig. G^b), and the 
pair correlation function g(r) has many peaks and a correlation length of r corr = (23 — 30) a. 

The particles were imaged through the top window by a video camera. We digitized 
movies of 256 frames at 30 frames per second. The 23.1 x 17.3 mm field of view included 
370 — 770 particles, Fig. [IJa). Coordinates x, y and velocities u x , u y were then calculated [f| 
for each particle in each frame. 

To apply a shear stress in a planar Couette configuration, in an ordinary fluid one uses 
moving plates, but in our dusty plasma, Fig.EJc), we used counter-propagating laser beams, 
as in Ref. Q], but with more laser power. Particles are pushed by the radiation pressure 
force F, which was verified [9] to be proportional to an incident laser intensity. We chose 
to use two opposing laser beams, rather than one, to avoid introducing a net momentum to 
the particle suspension. An Ar + laser beam, with a power that was varied up to 3.41 W, 
was split in two. At their foci at the particle location, the laser beams had a diameter of 
0.61 mm at the 1/e 2 level. Two rapidly oscillating scanning mirrors rastered the beams into 
vertical sheets, which made footprints on the lattice by striking it at an angle of 6.7 — 8.8° 
with respect to a horizontal plane. The sheets extended beyond the edges of the suspension. 

In discussing our results, we first note that steady-state conditions in the velocity profiles 
developed within 3 s after turning on the laser. All the maps and profiles of particle velocity 
and temperature presented here were computed well after that time. 

The particle velocity was always highest within the laser footprints. The velocity dimin- 
ished elsewhere, with a profile that we will model. 

Beginning with an undisturbed lattice, we applied increasing levels of shear stress, and 
we observed that the particle suspension passes through four stages elastic deformation, 
defect generation while in a solid state, onset of plastic flow, and fully developed shear 
flow. We present data for the latter two stages. At the onset of plastic flow, Fig. Hl[c), 
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the particles hopped between equilibrium lattice sites. Domain walls developed, and they 
moved continuously. The crystalline order of the lattice in the shearing region deteriorated, 
broadening the peaks in the static structure factor S(k) in Fig.^d). At still higher levels of 
shear stress, the lattice fully melted everywhere, and a shear flow developed, Fig. de). The 
particle motion was highly irregular on a small scale compared to the interparticle spacing, 
but on a larger scale, it was like a laminar flow in a fluid. The liquid-like order of the particle 
suspension is clearly indicated by the diffusiveness of the structure factor 5(k) in Fig. dT). 
Particles were confined so that after flowing out of the field of view they circulated around 
the suspension's perimeter and re-entered the field of view, and the suspension did not buckle 
in the vertical direction. Within the field of view, more than 95% of the time-averaged flow 
velocity was directed in the x-direction, with less than 5% of the flow velocity diverted in 
the ^/-direction. 

We used a method of binning the velocity measurements u x and u y according to a parti- 
cle's y position to compute spatial profiles v x (y) for the flow velocity, and the fluctuations 
parameterized as the kinetic temperature T x (y) and T y (y). Within each bin, we computed 
moments of the velocity, and then time-averaged over 256 consecutive frames, beginning 
after a steady state was attained. This yielded T x (y) and T y (y) from the second moments 
of u x and u y , respectively, as shown in Figs. Hfc) and[Tfe); and v x (y) from the first moment, 
as shown in Fig. Etc). We used the region between the laser sheets to calculate a spatially 
and temporally averaged < T y >, which we used to compute F y . We use F y and not F x 
to characterize our system, because there was no external energy input in the ^-direction. 
Separately, we also computed maps of velocity, u x (x, y), which were time-averaged to reduce 
noise, Figs. |2fa) and|U^b). 

At the onset of plastic flow, the velocity maps u x (x,y) reveal interesting 120° zigzag 
features, Fig.|2fa). This angle corresponds to a triangular lattice. The underlying movement 
of individual particles is best seen by viewing the movie jl^ . 

In the fully melted condition, the velocity map in Fig. [21(b) has no significant variation 
in the x direction. In this shear flow, particles slip past their neighbors 14]. 

The flow velocity profile v x (y) in Fig. |2fc) was curved, unlike in a traditional planar 
Couette flow where v x (y) is linear with y. We attribute this curvature to the frictional 
drag exerted on particles by the gas. A balance between this drag in the x direction and 
the viscous transport of particle momentum mu x in the y direction away from the laser 
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footprints accounts for the observed steady-state velocity profile. 

We modeled the velocity profile in the continuum approximation, using the Navier-Stokes 
equation with a term for the gas drag dw/dt + (vV)v = — p _1 Vp + (n/p)V 2 v + [(/p + 
?]/(3p)]V(V • v) — VdV. Here, the parameters for the continuum representing our particle 
suspension are v, p, p, 77, and (, which are the velocity, pressure, areal mass density, shear 
(dynamic) viscosity, and second viscosity, respectively, and v& is the gas friction. The ratio 
rj/p is the kinematic viscosity, which has the same dimensions in 2D and 3D systems. Our 
flow has a symmetry d/dx = and v y = 0. The Navier-Stokes equation is then reduced to 
d 2 v x (y)/dy 2 - {v d p/rj)v x (y) = and p = const. 

We will compare our results to a theoretical velocity profile v x (y) = [(Vi + V 2 e 2ah )e ay — 
(V2 + V1 e 2ah )e~ ay ]/ (e 3ah — e~ ah ), where a = Judp/r] and 2h is the distance between the laser 



sheets. We used the boundary conditions v x (—h) = —V\ and v x {h) = V 2 \ if V\ — V 2 — V, 
the theoretical profile simplifies to Vsmh(ay)/smh(ah). Here we have ignored the spatial 
dependence of rj/p, thereby neglecting the temperature gradient in the shear flow. As Ua — > 0, 
the solution approaches a linear velocity profile, as for a planar Couette flow. The gas drag 
Vd 7^ gives the steady-state velocity profiles a curvature depending on viscosity; this allows 
us to calculate the viscosity from the profiles, if v& is known. 

We fit our experimental velocity profiles in Fig. |2^c) to theory using experimental values 
of V\ and V 2 and a single free parameter a. We then calculated the kinematic viscosity 
t]j ' p— VdOT 2 using the known value of Epstein gas drag Ud = 0.87 s _1 for our experimental 
conditions 0. The resulting curves fit our profiles well. 

Our main experimental result is the kinematic viscosity rj/p of our particle suspension, 
Fig. OK a). Its value is of order 1 mm 2 s _1 , which is comparable to rj/p for both a 3D Yukawa 
system and liquid water |ll 0, S, 0, S| . For a given value of k, the parameter we varied was 
the applied laser power; increasing this power caused the shear stress to increase and r y to 
decrease, so that these two parameters were not varied independently. 

A prominent feature of Figs. Ota) and Otb) is a broad minimum in the viscosity. The 
minimum occurs in the range 70 < F y < 700. In 3D Yukawa systems, there is also a 
minimum, and the value of r for this minimum depends on k Q, O^y, Q]. 

We next consider some of the assumptions in our model and validate them using our 
experimental results. The weak dependence of viscosity on F y in the range 70 < F y < 700 
validates our neglecting the temperature gradient, but only in this range of T y \ our data 
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points outside this range are less reliable, due to this assumption. The Reynolds number of 
our shear flow was R = Vih/(rj/p) = 0.7 — 17; values this low typically indicate a laminar 
flow, validating that assumption in our model. 

We have also assumed that the neutral gas does not flow and entrain the particles. This 
is true because the gas-gas collision mean free path of 15 mm was of the order of the entire 
suspension's size, and the particles filled only a tiny solid angle. Thus, a gas atom that was 
struck by one moving particle was unlikely to strike another. 

The good fit in Fig. Efc) suggests that the Navier-Stokes model, which is a continuum 
model that does not describe motion of individual molecules, works well even when the ratio 
of the shearing region width 2h to the Wigner-Seitz radius a is as small as 17 to 24, as it 
was in our experiment. In Ref. |l|], tests suggested that validity of the Navier-Stokes model 
requires that the length scale 2h of the system exceed both the interparticle mean free path 
A m f p and r corr . Both of these conditions are satisfied for our experiment. We validate the 
2h > A m fp condition by observing that 2h/a was in the range 17 — 24 in our experiment, 
and judging that A m f p ~ a, based on the movies of the particle motion [l^j]. We validate the 
2h > r corr condition using the pair correlation results, yielding 2h/r corr = 3.0 — 9.4 for all 
data points except for two with 2h/r corr = 2.5 — 2.9; for these, F y > 1000. The latter two 
data points, which are for a mostly solid suspension, are therefore less reliable. 

However, the Navier-Stokes equation cannot be used inside the laser footprints, where 
the shear stress is applied, because this region is very narrow, less than an interparticle 
spacing. For this reason, we did not extend our fits into the footprints. In general, if the 
Navier-Stokes equation fails to apply, one cannot expect the usual relation, that the shear 
stress is the product of the dynamic viscosity and shear rate, to apply. 

Our experiment suggests a need for a theory or simulation for the viscosity of a 2D 
Yukawa system. The only previous work, to our knowledge, is for a 3D Yukawa system 

After we submitted this paper we learned of the experiment of Ref. [15j, where the shear 
viscosity of a multi-layered dusty plasma was measured using a shear flow induced by a 
single laser beam. 

We thank B. Liu, J. Marshall, R. Merlino, and F. Skiff for valuable discussions. This 
work was supported by NASA and DOE. 
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TABLE I: Parameters of the particle suspension measured without any shear flow, for three exper- 
iments. 

12 3 



Wigner-Seitz radius a (mm) 0.402 0.449 0.580 
k 0.36 0.42 0.53 

Q/e -11 940 -13 840 -16 360 

u P d = [Q 2 /2ire ma 3 ] 1 / 2 (s^ 1 ) 49.2 48.2 38.8 
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FIG. 1: Particle suspension at three levels of shear stress applied by a pair of counter-propagating 
laser sheets in a planar Couette configuration. Panels (a), (c), and (e) are top views with shear 
stress that is zero, moderate, and high, respectively. The snapshot of the undisturbed suspension 
(a) shows a highly ordered triangular lattice. At the onset of plastic flow (c), particle trajectories 
reveal localized hopping. In a fully developed shear flow (e), the trajectories show particles moving 
significantly everywhere. Profiles of the inverse particle temperature and T~ l are shown in the 
insets. The static structure factor S(k) was computed as the Fourier transform, with a Hanning 
window, of the raw bitmap images as in (a). The peaks in S(k) are distinctive for an undisturbed 
lattice (b); they are broadened and then diffuse at progressively higher levels of shear (d),(f). The 
particle suspension's parameters here and in Fig. |2] correspond to Experiment 1 in Table 1. We 
calculated rj, < T y >, and F y using data in the range —h<y<h. 



FIG. 2: Particle velocity data. Map (a), averaged over 64 frames at the onset of plastic flow, shows 
particle movement localized near the laser footprints. Map (b), averaged over 256 frames at the 
highest laser power, shows a fully developed shear flow that has no variation with x. Experimental 
profiles (data points) (c) are with fits to the Navier-Stokes model (curves), as used to calculate the 
kinematic viscosity of the particle suspension. 



FIG. 3: Results for kinematic viscosity rj/p, as functions of the shear stress (a) and coupling 
parameter F y (b). The latter two parameters were not varied independently, but were related 
as shown in (c). There is a broad minimum in viscosity in the range 70 < F y < 700. In (b), 
to aUow comparison to theory, we nor^ed the vtscosity as „• = „/(^>, where ^ , 
[Q 2 /27reoma 3 ] 1 / 2 for two dimensions [1 1| . 
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